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Abstract 



At the upper endpoint of the photon energy spectrum in T ^ X7, the 
standard NRQCD power counting breaks down and the OPE gives rise to 
color-octet structure functions. Furthermore, in this kinematic regime large 
Sudakov logarithms appear in the octet Wilson coefficients. The endpoint 
spectrum can be treated consistently within the framework of a recently de- 
veloped effective field theory of collinear and soft particles. Here we show 
that within this approach the octet structure functions arise naturally and 
that Sudakov logarithms can be summed using the renormalization group 
equations. We derive an expression for the resummed energy spectrum and, 
using a model lightcone structure function, investigate the phenomenological 
importance of the resummation. 
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I. INTRODUCTION 



Early theoretical analyses of heavy quarkonium decay were based on the color-singlet 
model (CSM). The underlying assumption of this model is that the heavy-quark-antiquark 
pair has the same quantum numbers as the quarkonium meson. (For example the bb that 
forms an T must be in a color-singlet ^Si configuration.) One consequence of such a restric- 
tive assumption is that theoretical predictions based on the CSM are simple, depending on 
only one nonperturbative parameter. The quantities first calculated in the CSM were the 
inclusive rates for quarkonium to decay into leptons and into light hadrons [0] . Subsequently 
the direct photon spectrum in inclusive radiative quarkonium decays was calculated 0. 

In recent years the simple CSM has been superseded by a nonrelativistic effective theory 
of QCD (NRQCD) [^Q. Inclusive decays of quarkonium are now understood in the frame- 
work of the operator product expansion (OPE), supplemented by the power-counting rules 
of NRQCD. In this formalism the direct photon spectrum of T decay is 

^ = 5:a(M,z)(T|0,|T), (1) 

i 

where z = 2E^/M, with M = 2mh. The are short-distance Wilson coefficients which can 
be calculated as a perturbative series in ^^(M), and the Oi are NRQCD operators. NRQCD 
power-counting rules assign a power of the relative velocity v of the heavy quarks to each 
operator and organizes the series. The series may be truncated at any order with omitted 
terms suppressed by powers of v. For S-wave mesons the formally leading-order contribution 
is the color-singlet operator, which is related to the wavefunction at the origin. Thus for 
5- wave decays one recovers the CSM at leading order in f . At higher orders in v color-octet 
operators need to be included. 

However, the picture of the photon spectrum in T — X7 decay which emerges is much 
richer than the naive expectation that the color-singlet contribution is leading 0-0]. At low 
values of the photon energy, fragmentation contributions to r(T X7) are important 
The situation at large values of the photon energy is even more interesting, because both 
the OPE and perturbative expansion break down. The breakdown of the OPE was first 
addressed in Ref. It was shown that the color-octet contributions, which give rise to a 
singular contribution at maximum photon energy, become leading for large photon energies. 
The singular nature is smeared by a nonperturbative structure function, which tames the 
endpoint behavior of the photon spectrum. The breakdown of the perturbative expansion 
gives rise to so-called Sudakov logarithms which have to be resummed. In a recent work 
it was pointed out that the leading Sudakov logarithms cancel in the CSM. However this is 
not the case for Sudakov logarithms in the color-octet contribution. 

Both the breakdown of the OPE and the appearance of Sudakov logarithms are symptoms 
of the same disease: NRQCD does not contain the correct low energy degrees of freedom 
to describe the endpoint of the photon spectrum. It does not contain collinear quarks and 
gluons. A theory constructed from the appropriate degrees of freedom was developed in 
Refs. [^|1T|. In those papers the theory was applied to the decay of a single heavy quark 



to light degrees of freedom. It was shown that the renormalization group equations (RGEs) 
in this theory sums Sudakov logarithms. In addition, for inclusive decays at the endpoint. 
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the nonperturbative structure function arises naturally from a modified version of the OPE. 
Here we apply the theory to the color-octet contributions to radiative T decay. In Section II 
we discuss the leading contributions in the endpoint region and motivate perturbative and 
nonperturbative resummation. In Section III we sum Sudakov logarithms using the RGEs 
in an effective field theory. In Section IV we introduce a phenomenological model for the 
shape function, convolute it with the resummed spectrum, and show how this changes the 
color-octet contribution to the spectrum. In Section V we conclude. 



II. LEADING ORDER RESULTS 

The inclusive radiative differential decay rate of T can be calculated using the optical 
theorem. This relates the decay rate to the imaginary part of the forward matrix element 
of the time ordered product of two currents 

f ^^L.(T|I,„r|T), (2) 

where we have used nonrelativistic normalization for the states: (T(P')|T(P)) = 
(2tt)^6^{P' — P). For large momentum transfer the time ordered product can be expanded 
in terms of local operators giving 

^ = j:a{z){T\0,\T). (3) 

The Wilson coefficients, Ci{z), can be calculated perturbatively as a series in as{M). The 
parametric size of the long distance matrix elements are determined by power counting in 
NRQCD, but to obtain quantitative results these matrix elements have to be extracted from 
experiments or lattice calculations. 

At leading order in the NRQCD v expansion, only the color-singlet, spin-triplet operator 

O^CS,) = E[<^^X-p'] [xtp^Yp] (4) 
p,p' 

contributes, with the leading order Wilson coefficient 

(n) X , . 128 a Co 



\2 



(5) 



where cq = -1/3 for T. This is the CSM result 0, with {T\Oi{^Si)\T) (3/27r)|i?(0)|2, 
where -R(O) is the radial wavefunction at the origin. The first color-octet contributions to 
(I) are suppressed by f ^ relative to this color-singlet contribution. There are two operators 

p.p' 

OsCPo) = ^ E P • ^ T^X^P'] [X^ P ■ o- T>p] , (6) 
p>p' 



3 



with leading order Wilson coefficients ^ 

cf(z)=C'f5(l-^) 



(7) 



where 



16 as a egTT 

448 as a CqTi 
M4 



(8) 



Since the color-octet Wilson coefficients are enhanced by a power of as{M)/TT relative to the 
color-singlet one, the overall suppression of the color-octet contribution is t^v'^ / as{M) ~ w^, 
where we have used that numerically as/^r ~ f^. 

The singular nature of the coefficients (|^) is an indication that the OPE is breaking down. 
We can obtain a rough estimate for the value of z at which the octet contributions become 
of order the color-singlet contribution by smearing the perturbative spectrum over a small 
region near z = 1. Integrating over 1 — f ^ < 2; < 1 gives a color-singlet contribution that 
scales as a1{M)v'^ and a color-octet contribution that scales as as{M)T[v'^. Thus, the ratio 
of octet to singlet in this region of phase space is 7rt>^/as(M) ~ 1, making the color-octet 
contribution of the same order as the color-singlet one. 

It was shown in Ref. that in precisely this endpoint region the OPE breaks down 
and an infinite set of operators have to be resummed into lightcone distribution functions 
Each structure function gives the probabihty to find a bb pair with the appropriate 
quantum numbers and residual momentum fc+ in the T. For the color-singlet contribution 
the structure function can be calculated using the vacuum saturation approximation. It 
simply shifts the maximal photon energy from 2m;, to Mf . The color-octet contributions, 
however, give rise to two new nonperturbative functions at leading order. At higher order 
there are an infinite number of additional structure functions, so the differential decay rate 
in the endpoint region is 



dV 

dz 



ldk^Qiz,k+)f,ik+){T\0,\T) 



(9) 



Another effect one encounters in the endpoint region is the appearance of Sudakov log- 
arithms in the Wilson coefficients, which ruin the perturbative expansion. Consider, for 
example, the Wilson coefficients for the color-octet operators at next-to-leading order in as- 
In the z ^ 1 limit they are 



C^s\'So)iz) = ^C^^\'So) 




23 _nf_ 
6^3 



1 



1 - z 
1 

1 - z 



(10) 



^This is very similar to the behavior of the OPE for B — > Xuiv at the endpoint of the lepton 
energy spectrum and B Xg^ at the endpoint of the photon energy spectrum |12|. 
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If these coefficients are integrated over the shape-function region (1 — < 2 < 1), then the 
first plus distribution on the right-hand side gives rise to a double logarithm, log^f^, while 
the second plus distribution gives a single logarithm, logf^. Both of these are numerically 
of order l/a^. This clearly ruins the perturbative expansion. Therefore, to obtain a well 
controlled expansion, these logarithms must be summed. 

III. SUMMING SUDAKOV LOGARITHMS 

The NRQCD power-counting rules break down as z — 1 because NRQCD does not in- 
clude all of the long distance modes: coUinear physics is missing from the theory. An effective 
theory which includes coUinear physics was developed in Refs. [0,0. This coUinear-soft 
theory describes the interactions of highly energetic coUinear modes with soft degrees of 
freedom. To describe T decay at the endpoint we have to couple the coUinear-soft theory 
with NRQCD This is analogous to 5 — > decays at the endpoint, which was studied 
in the context of effective field theory in Ref. 0]. We will closely follow the development 
in that paper. 

Understanding inclusive T X7 decays near the endpoint is a two-step process.0 In 
the first step we must integrate out the large scale, M = 2mh, set by the bb pair constituting 
the T. This is done by matching onto the coUinear-soft theory. In the second step coUinear 
modes are integrated out at a scale which is set by the invariant mass of the coUinear jet. This 
is done by performing an OPE and matching onto a soft theory containing operators which 
are nonlocal along the lightcone and whose matrix elements are the lightcone distribution 
functions, fi{k^). Sudakov logarithms are summed by using the effective theory RGEs. 
Operators are run from the hard scale M to the coUinear scale where the OPE is performed. 
The nonlocal soft operators are then run from the coUinear scale down to a soft scale where 
their matrix elements do not contain large logarithms. This procedure sums all Sudakov 
logarithms into the short distance coefficient functions. 

To better understand the scales involved consider the momentum of a coUinear particle 
moving near the lightcone. In lightcone coordinates we can write this momentum as p = 
{p'^ ,p~ ,p±)- Since the mass of the particle is much smaller than its energy, we define 
p2 M^A^, where M is the scale that sets the energy and A is a small parameter. The 
lightcone momentum components are widely separated. If we choose p~ to be 0{M), then 
P_l/p~ ~ A, and p^ /p" ~ A^. We refer to these two scales as coUinear and soft, respectively. 
To be concrete consider the bb pair to have momentum Mf ^ -|- k^^, where f ^ = (1, 0, 0, 0) and 

is 0{Aqcb) in the T center-of-mass frame. The photon momentum is Mzn^/2, where 
we have chosen n'^ = (1,0,0,1). In the endpoint region the hadronic jet recoiling against 
the photon moves in the opposite lightcone direction = (1,0,0,-1), with momentum 
Px = Mn^/2 + M(l — z)n^/2 + k'^ . Thus the hadronic jet has n ■ px = Px ^ Next note 
that m\ ~ ^^(1 - z). For {I - z) ^ ^ Aqcd/M we find 



2t 



In Ref. [13 1 it was pointed out that at higher orders in perturbation theory one has to adopt a 
one-step scheme similar to the one developed in a slightly different context in [0|. 
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FIG. 1. Matching onto operators in the effective field theory with one and two gluons in the 
final state. 



mx ~ VMAqcd , (11) 
which is the colhnear scale. This imphes that for this process the colhnear-soft expansion 



parameter A is of order a/1 — z ~ ^Aqcd/M. The soft scale is the component of the 
hadronic momentum in the n direction: 

2 

n-px ^ ~ Aqcd ~ MX^ . (12) 

n-px 

Thus in order to sum Sudakov logarithms in T — > X7 we first match onto the coUinear-soft 
theory at M and run operators in this theory down to the collinear scale /ic ~ Ma/1 — z. At 
that scale we perform the OPE by matching onto a soft theory containing operators that are 
nonlocal along the lightcone and run these operators down to the soft scale /is ~ M(l — z). 

A. The collinear-soft theory 

We first need to integrate out the large scale M by matching onto the collinear-soft 
theory. This is done by calculating matrix elements in QCD, expanding them in powers of 
A, and matching onto products of Wilson coefficients and operators in the effective theory. 
For the process of interest this matching is illustrated in Fig. |I[ The heavy-quark spinors and 
the heavy-quark propagator in QCD can be expanded in powers of v to match onto NRQCD. 
The QCD spinor is decomposed into two two-component Pauli spinors if) and x for th^ heavy 
quark and anti-quark, respectively [0.0 We also need to expand the amplitude in powers 
of A to match onto the collinear-soft theory. This is done by using the power-counting rules 
for the gluon field given in Ref. [|lT| and by scaling the different components of the gluon 
momentum as 



0(1) + 0(A) + 0(A2). (13) 



Omitting the straightforward but unenlightening details, the color-octet contributions match 
in the T rest frame, at leading order in A, onto the operators 



Note that in Ref. |14] states are normalized relativistically. 
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FIG. 2. Collinear and soft diagrams needed to calculate the renormalization of the vertices in 
the effective theory. 



p 

Q^sCPj) = Ag,e, [gl'g'^'' + gT g'' - gTn^n') A ■ P.[xlpA • ^sT^^] , (14) 

p 

where e"^ = e°'^''^npVp and g^ = g^^ — {n^n'^ + n'^n^)/2. Hatted variables are divided by 
M, and the A'^* are boost matrices which boost from the hh center-of-mass frame [Q. We 
have chosen factors in the operators such that at leading order the corresponding Wilson 
coefficients satisfy 

C^{^S,) = C^{'P,)^C^ = 1. (15) 

Note that there is no color-singlet operator at this order in A, and therefore leading Sudakov 
logarithms are absent in the CSM 

To calculate the renormalization group equations of these operators, we also need the 
Feynman rules shown in Table |I[ They are all obtained by expanding full theory diagrams in 



powers of A and v. In addition the coupling of soft gluons is given by HQET and LEET |T5 



Feynman rules, and the coupling of three collinear gluons is identical to the three gluon 
vertex in QCD. The renormalization is independent of the operator and the two diagrams 
that need to be calculated are shown in Fig. H. The result for the collinear and soft graphs 
are 




A = (ZT^iy^) ^ + + 2e)r(l - 2e)i Vr , (17) 

where K"''" = iy^^\^So),V^^\''Pj)). 

To renormalize the vertex, we expand in e, keeping only the divergent pieces. This 
must equal {Z — ViV,""^"" . Note however that Z is not the counterterm for the vertex, 
rather Zq = Zyj^Z^"^ Z~^ Zl^l''^ ^ where Z^ = \ since we are not considering QED corrections, 
Zy^j^ is the counterterm of the color-octet hji^ current, and Zj, is the gluon wave function 
counterterm: 



1 5 2\ 
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TABLE I. Feynman rules in the collinear-soft theory at leading order in A. The vertices 
^i,8('^5'i), Vi^8(^-Pi) have a zero matching coefficient at this order. 

Diagram Feynman Rule 



k, [i, a 



p, \, b 



i 




q, p, c 



1^ jafec^M(2^'^Pn ■ p - n^pP - nPp'' - n^nPfi ■ k) 




^8 y ^oj 



a, a 




P2, (3, 



- 1^ [gf'lg^'^ + gf'lg^' - g^^rf, 
xA-p, r^lpA • o-^T^^p 



This leads to 



4:71 



17' 



M2 



6e 



3e 



(19) 



We can check this result by matching the effective theory to QCD while regulating all 
divergences using dimensional regularization. In this approach there is no scale in effective 
theory loop integrals so they are zero. This leaves only the counterterm Z, which must 
match the e poles in the QCD calculation. The QCD vertex at one loop can be extracted 
from a calculation by Maltoni, Mangano, and Petrelli [|T6|. Equating the pole terms in the 
QCD calculation to Z we again obtain (|19D. 

The RGE for the Wilson coefficients of the operators in (|1^ is 



(20) 
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FIG. 3. The leading OPE: tree level matching of the time ordered product in the collinear-soft 
theory to a nonlocal operator in the soft theory. 



in the same form as in Ref. [p7 



In order to make use of previous results from B Xg'y, we write the anomalous dimension 

(21) 



rrih 



Defining 



cusp 



TT 



TT 



-Bi, 7 = 7i, 



TT 



TT 



we find 



rf = Ca, B, + 271 



-Ly A ^ 



(22) 



(23) 



where /3o = (IIC^ — 2nj)/3. At this point we can only determine the above linear combi- 
nation of Bi and 71, and we do not know the expression for r2'^-'. However, as we will show 
below, Bi and T'^'^^ can be found in the soft theory. To sum the leading Sudakov logarithms 
in the collinear-soft theory, we run the operator given in ([l^) from the matching scale /i = M 
down to a scale /ic ~ wMAqcd at which the OPE is performed. We can lift the solution 



from Ref. |TT 



log 



Cv{M) 



adj 



/3o'«s(M) 
Bi + 271 



1 + log ?/ 



PI 



I ~ y + y\ogy - - log^ y 



logy - 



PI 



y-l-\ogy 



(24) 



where y = as{ij) / as{M). 



B. The purely Soft Theory 

At the collinear scale /ic ~ Ma/1 — z we integrate out collinear modes and perform an 
OPE for the inclusive T radiative decay rate in the endpoint region. The result is a nonlocal 
OPE in which the two currents are separated along a light-like direction. Diagrammatically 
this is illustrated in Fig. |^. We write the momentum of the jet as 

Px = + + ^(1 - x)n^^ , (25) 
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where k'^ is the residual momentum of the bb pair. Note we distinguish x from z, because 
the momentum of the jet is not exactly the same as the momentum of the collinear gluon 
which was integrated out. The two can differ slightly due to the emission of soft quanta 
by the jet. The imaginary part of the tree level diagram on the left hand side of Fig. ^ is 
proportional to S{p1). Taking k'^ ~ M(l — a;) ~ MA^ and expanding in the small parameter 
A, we match at leading order onto the operator 

0^ix) = EK'Hx-p'] 5(1- x + ib+) [xLpFi^p] , (26) 
pp' 

where r8(^S'o) = and r8('^Po) = T'^p-cr/v^. The x serves as a continuous label on the 
operator. 

Each operator has a Wilson coefficient, which can depend on the photon energy fraction z. 
The differential decay rate is given by a convolution of the matrix elements of the operators 
Oi{x) and the corresponding Wilson coefficients 

— = Y^jdx C,{x - z- fi)Mx- /x)(T|0,|T) , (27) 

i 

where 

are the lightcone distribution functions. 

The convolution of the short distance Wilson coefficients and the long distance operators 
presents a technical problem since the RGE for 0{x) will be given in terms of a convolution as 
well. We use a Mellin transform to deconvolute (P7D, which is equivalent to taking moments 
of the decay rate. We must restrict ourselves to large moments A^, since it is the limit 
N ^ oo which is equivalent to the region z ^ 1. Once the the final expression in moment 
space is obtained we can take an inverse-Mellin transform to get back to 2;-space. This 
procedure is valid up to corrections of order 1 — z. Taking large moments of the expression 
in (P7| ) gives 

r dV 

T{N) ^ j dzz''-'^=J2QiN;f^)f,iN;f^), (29) 

where 

MN-fi)= /a;^-V,(x;/i). (30) 



To match onto the soft theory, we compare large moments of the differential decay rate 
calculated in the coUinear-soft effective theory and the soft effective theory in the parton 
model. Large moments of the one loop expression calculated in the collinear-soft theory are 
given in ( |A7| ). The one loop expression for {bb\0{N] fi)\bb) can be lifted from Ref. [|TU| with 
the replacement — > Ca- 



{bb\0{N;f^)\b-b) = l-^ 



2 /"iV fiN 

41og — 41og— — 

Mrio Muq 



(31) 
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where no = e At the scale /Xc = M^tiq/N all logarithms match, and at that scale the 
tree level matching coefficients are 



Ci{N■^^) = C^ 



(0) 



(32) 



In the matrix element (|3T|) all logarithms vanish at the scale /i^ = Muq/N. To sum the 
large logarithms, we therefore have to run the Wilson coefficient in the soft theory, Ci{N] fi), 



from /Xc = M^tiq/N to /i^ = Muq/N. Again, we keep the notation introduced for B Xg'y 
and write the RGE as 



^-^CiiN;fi) = ^{N;fi)CiiN;fi), 



with 



71 



iN;fi) = 2r:tiM\og^ + B 



MriQ 



(33) 



(34) 



From the results in Refs. |T0|,|T8|,|T9| and the previous section we obtain 

'67 7r^\ 



= Ca, 
Bi = —Ca 



padj 
-L 2 



Ca 



Ca 



36 12 



(35) 



where we have taken from Refs. IIT8|,p!9 



The solution to (^3|) combined with the running in the collinear-soft theory (|^) can be 
lifted directly from Ref. by substituting into that result the expressions in (pSl). This 
gives the fully resummed result in moment space 



V{N) = Y.Cf'^f{N-Mno/N)e 



log{Ar)gi(x)+92(x) 



(36) 



where 

9i{x) 
92{x) 



2Tt 



adj 



2_ 

Pi 



[(1 - 2x) log(l - 2x) - 2(1 - x) log(l - x)] 



g-padj 

' ^ -log(l-2x)+21og(l-x)] 



2rfj/3i 



Pi 



log(l - 2x) - 21og(l - X) + 2 log'(l - 2X) - log'(l - x) 



471 



25i 



4ri 



adj 



log(l-x) + -7^1og(l-2x)- ^ 

Po Po Po 



log no [log(l - 2x) - log(l - x)] , (37) 



X = log(iV) a,(M)/3o/47r, and /3i = (34^2 _ lOC^nj - 6C^n/)/3. 
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IV. RESULTS 



Now that we have the resummed rate in moment space (P^D, we must take the inverse- 
Melhn transform to obtain the expression for the photon energy spectrum. Fortunately, the 
inverse-Melhn transform of the resummed rate in B ^ Xg'-f decay has been calculated in 
Ref. and we can use that result by simply substituting in the gi from (0). We obtain 
the following for the resummed Wilson coefficients of the octet contribution 

r(o) ^ fzi/^ .^ exp[/^i[a,^o//(47r)] + c/2MoV(47r)]] | 

^ ' ^4 ^ ^ r[l-g^[a,Pol/W]-asPol/Wg[Ml/W]]j ' ^^^^ 

where (7^°^ = Cg\^So), Cg\^Po), and I ^ — log(x — z). Each Ci{x — z) is evaluated at 
the soft scale so that all leading and next-to-leading Sudakov logarithms have been summed 
into it. 

One way of checking (^) is by expanding in powers of and comparing to the fixed 
order calculation. Using 




6(1- z- r,)-^^ — ^ + 5(1 - z) ^ 

^ " l-z ^ ' n+l 



(39) 



as the definition for plus distributions it is straightforward to verify that the order as term 
in the expansion of (|38D reproduces the plus distributions in (|T0|). 

Recall that the covariant derivative in ( |26| ) is of order Kqqd/M. If we consider the limit 
1 — z ^ Aqcd/^; then the covariant derivative in the operator appearing in the structure 
function can be neglected. In this limit we can perform the integral in (^) to obtain 

l-EC.d-.). (40) 

This result gives the effect of the perturbative resummation without the structure function. 
The quantity Ci{l — z)/Cf'\ which is the same for the two leading octet configurations, is 
shown as the dashed line in Fig. ^.0 

However, for the bb system mv"^ ~ Aqcd so the covariant derivative cannot be dropped 
in the endpoint region of the photon spectrum, and the differential rate is given by the 
convolution in (^). The lightcone distribution function is a nonperturbative function and 
needs to be modelled. In this paper we will be content with the simple structure functions 
introduced in Ref. for inclusive B decays 

/(fc+) = iV 1^1 - e(^+'^)'=-'/^ , (41) 

where is chosen so that the integral of the structure function is normalized to one. In 
principle the structure function can be different for the different color-octet states. But 



^The Landau pole in ( |38D should be dealt with in the same fashion as in the B decays pl] |. 
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z 

FIG. 4. The differential decay spectra near the endpoint region 0.7 < z in arbitrary units. The 
solid curve is the perturbative resummation convoluted with the structure function and the dashed 
curve is the perturbative resummation Cj(l — z)/cf\ The dotted curve is the plus distribution 
terms in the one-loop result JT^j, and the dot-dashed curve is these terms convoluted with the 
structure function. 



since we are ignorant of the nonperturbative structure function, we will naively use the 
same model for both the ^5*0 and ^Pq configurations. The structure function for B mesons 
have the property that the first moment with respect to vanishes. For quarkonium the 
first moment of (|28D (where x = 1 + k^) with respect to is 

^ ^ (T|Ep.p4xlp^T"F.^p,]^D+[^tT"F.X-p]|T) 
(T I Ep.p' [xlp'T^Fi^p,] [^T'^F.X-p] I T) 



Therefore, we need to shift k'^ in pT| ) to k'^ + Ai, so that ( |41| ) will have the desired first 
moment for quarkonium decays. The integration limits for k~^ are, similar to the case for 
B decays, from —M to My — M. Both A and Ai are nonperturbative parameters related 
through A = Mf — M — Ki. We use the following numbers in our plots: = 0.2, a = 1/137, 
rub = 4.8 GeV, Mr = 9.46 GeV, a = 1, A = 480 MeV, and Ai = -620 MeV. 

In Fig. 1^ the convolution of Ci{x — z)/Cf'' with the model of the structure function (^T]) 
is shown as the solid line. In addition we show as the dotted line the terms in the NLO QCD 
expression that dominate in the endpoint region (p!0| ) divided by Cf'' , and as the dot-dashed 
line the convolution of these terms with (0). Thus Fig. |^ gives a picture of the effects of 
resummation. The singular plus distribution piece of the NLO QCD expression is tamed by 
both the perturbative resummation and the structure function. Either of these alone gives 
a similar curve, which is peaked near z = 0.94. However, to be consistent the perturbative 
resummation must be convoluted with the structure function. This gives a curve that is 
broader, with a peak that is 34% lower and shifted to 2; = 0.87. Changing the values of the 
structure function parameters changes the shape of the curve. Halving the value of A gives 
a narrower peak that is 30% lower and shifted to 2; = 0.83. If is increased by 10% the 
peak moves slightly to the left and decreases in height by 5%. Doubling the value of a in 
(PI) slightly raises the peak, and steepens the curve as it goes to zero at the endpoint. 
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FIG. 5. The differential decay spectra in the region 0.5 < z. The dashed curves are the fully 
resummed result convoluted with the shape function for two choices of the octet matrix elements. 
In addition we have interpolated the fully resummed result with the NLO result in the region away 
from the endpoint. The dotted curves are the NLO result convoluted with the structure function for 
two choices of the matrix elements. The solid curve is the color-singlet contribution. 



In Fig. 1^, we show as dashed curves the fully resummed color-octet contribution convo- 
luted with the structure function for two choices of the matrix elements. We also show as the 

(3^i)|T) = 3.63 GeV^ 



23 



solid curve the color-singlet contribution. Here we used (T|Oi(^S'i)|T) 
The values of the color-octet matrix elements are not well determined. One may be tempted 
to use a naive power-counting argument which gives (T|08(^>S'o)|T) ~ {T\Os(c'Po)\T) /ml ~ 
f^(T|(9i('^S'i)|T), with ~ 0.1. However, this gives values for the color-octet matrix el- 
ements that are too large to be compatible with data on the ratio of hadronic to leptonic 
decays 0. The data suggest that the octet matrix elements are at least an order of mag- 
nitude smaller than the naive power-counting estimate. Therefore, we use values that are 
10 times and 100 times smaller than estimates from naive power counting, with the larger 
matrix element yielding the higher peak.Q We have added to the resummed result the NLO 
QCD result with the singular terms (0) subtracted off. This interpolates between the NLO 
QCD expression at lower values of z and the resummed result in the endpoint region. For 
comparison we show as the dotted curves the NLO QCD contribution convoluted with 
the shape function for the two choices of the color-octet matrix elements. 

Since the octet contributions dominate the color-singlet contribution in this endpoint 
region (or at least are of the same order of magnitude) it should be possible to use this process 
to constrain the size of the color octet matrix element. The suppression of the matrix element 
compared to the naive power counting estimate makes a measurement of this matrix element 
particularly interesting and might shed some light on the convergence of the v expansion 



^In Ref. 1 24 1 it was argued that a factor of l/2Nc should be included in a naive estimate of the 
color-octet matrix elements. The larger of our two choices for this matrix elements is of the same 
order of magnitude as this modified naive estimate. 
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in NRQCD. However, before a meaningful comparison of theory to data on radiative T 
decay can be made, subleading Sudakov logarithms in the color-singlet contribution must 
be summed. The existence of these logarithms was first pointed out in Ref. where it was 
observed that though leading Sudakov logarithms cancel in the color-singlet contribution to 
the differential rate, they are present in the derivative of the rate at the endpoint. We leave 
the resummation of these logarithms to a future publication p5 . 



V. CONCLUSION 



Using an effective field theory approach, we have resummed Sudakov logarithms in the 
leading color-octet contributions to the T —>■ Xj differential decay rate in the endpoint 
region. This is done in two steps. First we match onto an effective theory with coUinear 
and soft degrees of freedom and run the theory to the coUinear scale. Next we integrate out 
collinear modes by performing an OPE, matching onto non-local operators which are run to 
the soft scale. This sums all Sudakov logarithms into Wilson coefficients of these operators. 
The color-octet contribution to the differential decay rate in the endpoint region is given by 
the convolution of the Wilson coefficients with matrix elements of the operators between T 
states. The latter are the color-octet structure functions defined in Ref. JTJ. 

We choose a simple model for the structure functions to investigate the phenomenological 
consequences of resummation. We find that either the perturbative resummation or the 
inclusion of a structure function cures the singular behavior of the QCD results. Both give a 
similar effect, causing the curve to turn over near z = 0.94 and to go to zero at the endpoint. 
However, the effective field theory approach makes it clear that the correct expression for the 
differential rate near the endpoint is given by a convolution of the perturbative resummation 
and the structure function. This gives a spectrum that has a broader and lower peak than 
we obtain by including only perturbative resummation or the structure function, shifting 
the peak to z = 0.87. 

Before a meaningful comparison to data can be made the color-singlet result must be 
resummed as well. Only then can data on the decay spectrum be used to constrain the size 
of the octet contribution in the endpoint region [ p5[] . 
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APPENDIX A: FORWARD SCATTERING MATRIX ELEMENT 



Matching the forward scattering matrix element in the QCD and coUinear-soft theory 
gives an important check that we are reproducing the infrared physics of the full theory. We 
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have already calculated the vertex corrections. Note that there are two contributions from 
the vertex loops so the vertex contribution to the forward scattering matrix element is twice 
that given in (|16|). There is a contribution to the forward scattering matrix element from a 
ladder graph, which we have not evaluated yet. It gives 



2n 



Att- 



r(l + e)r(l + 2e)r(l-2e)^A^o, 



where 



(0) 



pl-i5 



(Al) 



(A2) 



is the tree level amplitude, and = M^(l — z). In addition there is a contribution coming 
from virtual corrections to the coUinear gluon propagator, which include the fermion, gluon, 
and ghost loops. 
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Adding all contributions gives the full one-loop corrections to the forward scattering ampli- 
tude 
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- C^log^Kl - z)i-l - t6)] -Ca- log[(l -z){-l- i5)] 
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+ ^log[(l-^)(-l-25)]}A^o + ... 



(A4) 



where the ellipses represent terms of order not enhanced by logarithms. 

Next add the counterterms. There is a gluon propagator insertion that gives —M.q{Z^ — 
1). The vertex insertion is Zyj^Z^j'^ Z'q' . In order to facilitate the comparison with the 
results in Ref. 0, we use dimensional regularization to regulate the infrared divergences in 
the heavy quark sector, and therefore Zy^y^ = 1. Adding all the counterterms cancels the e 
poles in ( |A4|) , leaving 
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n 



Ca^ log[(l - z)i-l - t6)] + ^ log[(l -z)i-l- t5)] \Mo + .... 

using 



6 " 3 

To obtain the differential decay rate we take — l/vr Im 
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Comparing to Ref. (taking the limit z — > 1 in their expressions), we confirm that at 
the matching scale fi = M the effective theory reproduces the plus distributions in the full 
theory. Taking large moments of the imaginary part of ( |A5| ) gives: 
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